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Abstract. Using the formalism of extended = 4 supersymmetric quantum mechanics 
we consider the procedure of the construction of multi-well potentials. We demonstrate the 
form-invariance of Hamiltonians entering the supermultiplet, using the presented relation 
for integrals, which contain fundamental solutions. The possibility of partial A^ = 4 super- 
symmetry breaking is determined. We also obtain exact forms of multi-well potentials, both 
symmetric and asymmetric, using the Hamiltonian of harmonic oscillator as initial. The 
modification of the shape of potentials due to variation of parameters is also discussed, as 
well as application of the obtained results to the study of tunneling processes. We consider 
the case of exact, as well as partially broken A^ = 4 supersymmetry. The distinctive feature 
of obtained probability densities and potentials is a parametric freedom, which allows to 
substantially modify their shape. We obtain the expressions for probability densities under 
the generalization of the Ornstein-Uhlenbeck process. 

Key words: supersymmetry; solvability; partial breaking of A^ = 4 supersymmetry; stochas- 
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1 Introduction 

Hamiltonians for systems with multi-well potentials are in the focus of classical and quantum 
dynamics. One of the most fascinating applications of multi-well potentials in classical dynamics 
is the description of particles transitions from one local minimum to another under the influence 
of different types of noise. This problem is general for many physical systems from Universe 
to microphysics. Calculations of the transition rates from one local minimum to another were 
initiated by Kramers [1] and after 70 years of development the problem is far from complete [2] . 
From the modern perspective another effect of amplification of transitions over the barrier under 
the action of weak time dependent periodic signal, which is called stochastic resonance [3l [4] 
attracts the special attention. 

The quantum mechanical dynamics in multi-well potentials is interesting for its own reasons. 
We just recall the tunneling effect, which is in focus again in respect to the wide-ranging research 
of trapping of atoms of alkali metals in superfluid state [SI E]. Moreover, when more than 
one barrier exists, the resonant amplification of tunneling rate is possible, which is known as 
resonant tunneling [71 [5]. Experimental observation of this phenomenon in superconductive 
heterostructures forms a basis for construction of resonant tunneling diode. We have to point 
out the importance of having a mechanism of changing the parameters of the potential (locations 
of minima, heights of barriers etc.), to determine the conditions of resonant tunneling. 



*This paper is a contribution to the Proceedings of the Workshop "Supersymmetric Quantum Me- 
chanics and Spectral Design" (July 18-30, 2010, Bcnasquc, Spain). The full collection is available at 
|http://www.emis.de/journals/SlGMA/SUSYQM2010.html| 
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Classical and quantum systems mentioned in the above, are common in methods of their 
analysis. The dynamics of stochastic systems is described by the Fokker-Planck equation (FP) 
[9l[T0]. One of the methods for solving the FP equation is the eigenfunction expansion method. 
Due to a formal similarity of the FP and the Schrodinger equations, the eigenfunction expansion 
method, which is analogous to the bound state expansion of the Schrodinger equation, can be 
applied to the Fokker-Planck equation. For instance, in the case of bistable stochastic system, 
the knowledge of the full set of wave functions and eigenvalues of the corresponding quantum 
Hamiltonian completely determines the time evolution of the solutions to the FP equation. 
Moreover, it allows one to make conclusions about the dynamics of the corresponding metastable 
system. Analysis of processes in multi-well potentials is complicated due to the fact that the 
existing models deal usually with piecewise potentials (such as constructed from rectangular or 
parabolic wells and barriers), which are presumably far from the real potentials. That is why 
analytic wave functions and spectrum in such potentials are unknown, which implies the only 
numerical analysis of their properties. 

Solving the FP equation exactly in this approach is closely related to the existence of exactly 
solvable quantum-mechanic problems, the number of which significantly increases during the 
last years [111 [T2] . Recall that the first meaningful example of construction of the significantly 
nonlinear models of diffusion in bistable system was constructed in [13^ I14j. The formalism of 
the Darboux transformation, used in these papers, is the one of the basic methods in construct- 
ing the isospectral Hamiltonians in supersymmetric quantum mechanics (SUSY QM). Further 
development in the construction of the exactly solvable stochastic models was achieved in [15] , 
where the Fokker-Planck models with prescribed properties were constructed by use of the 
Darboux-Grum procedure. Existence of exactly solvable |16] and partially solvable [TTj quan- 
tum mechanical models with multi-well potentials could fill the gap in theoretical analysis of 
the above mentioned processes and may be considered as a more realistic approximation in their 
research. Approach based on the accounting the instanton contribution in double-well-like po- 
tentials \18\ [T9] is commonly used in the description of tunneling processes and, in particular, 
in studying the features of the Bose-Einstein condencates in multi-well traps. The instanton 
calculus is also well applied to the transition between the wells under the influence of noise, 
because the dynamics of these processes is mainly determined by the energy of the first excited 
state. 

This paper is aimed at determining the relations between different types of potentials in the 
framework of the extended supersymmetric quantum mechanics {N = 4 SUSY QM) [201 [211 [22] . 
We also study the possibility of changing parameters of potentials in wide range and obtaining 
the exact expressions for spectrum and wave functions. The latter is especially important when 
potential is equipped with two and even more local minima, so the "resonant" tunneling may 
be realized. The obtained expressions are used for derivation of new exactly solvable stochastic 
models. Here we develop the approach of j23j to the construction of exactly solvable stochastic 
models for potentials with several local minima. Contrary to the case of previously established 
connection between the FP equation and N = 2 SUSY QM, the considered approach allows one 
to extend the range of exactly solvable stochastic models, since the super Hamiltonian of = 4 
SUSY QM contains large number of the isospectral Hamiltonians. A special characteristic of 
the approach is the existence of a parametric freedom in probability densities and new poten- 
tials entering the Langevin equation. It makes, in particular, possible to change the shape of 
potentials and densities without changing the time dependence of the probability density. This 
fact makes possible to study the characteristics of stochastic dynamics of particles in multi-well 
potentials with variable parameters. On the other hand, the development of new exactly solvable 
diffusional models leads to the construction of more adequate approximations to real processes. 

In Section [2] we briefly discuss the procedure of construction of isospectral Hamiltonians with 
additional states below the ground state of the initial Hamiltonian in A^ = 4 SUSY QM. We es- 
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pecially emphasize that multi-well potentials can arise due to general properties of the solutions 
to the Schrodinger type auxiliary equation. In particular, we discuss the phenomenon of partial 
supersymmetry breaking in = 4 SUSY QM. In the next section we give general expressions, 
that allows one to analyze the obtained potentials and wave functions, independently on the 
concrete form of the initial Hamiltonian. It allows, for instance, to generalize the concept of 
form-invariant potentials [24] and to calculate the normalization constants of zero-modes wave 
functions in the case of exact and partially broken supersymmetry. Section U] is devoted to 
the construction of stochastic models, by use of connections in the supermultiplet of isospectral 
Hamiltonians of = 4 SUSY QM. We obtain the expressions for corresponding probability 
functions and potentials, entering the Langevin equation. Using the harmonic oscillator Hamil- 
tonian as the initial one, we obtain the expressions of the potentials of isospectral Hamiltonians 
and wave functions for the case of exact and partially broken supersymmetry. The analysis of 
conditions for emergence of double and triple-well potentials and ways for varying of their form 
in wide range is provided in Section [5j The proposed scheme of construction of new stochastic 
models is demonstrated in Section [6] on the example of generalization of well known Ornstein- 
Uhlenbeck process. We give the results of the calculations of probability densities and discuss 
the ways of their modification when parameters of scheme vary. In conclusion we summarize 
the main results and their possible applications. 



Isospectral Hamiltonians of AT = 4 SUSY QM 
with additional states 



Extended = 4 SUSY QM [20l [HI [22] is equivalent to the second-order polynomial SUSY QM 
(reducible case) [25l [26l [27l [28] and assumes the existence of complex operators of sypersym- 
metries Qi (Qi) and Q2 {Q2), through which the Hamiltonians H^^ can be expressed. Hamil- 
tonian of A'^ = 4 SUSY QM has a form (h = m = 1): 

^ ' .2 , T/2/ N , ^(1)t/'/^^A = 1/'^2 , T.2/^N , ^(2), 



V,{x) = W'{x) + lai 



1 (,) W"{x) 
2^3 w'{x) 



where W{x) is a superpotential and (jg*^ are matrices, which commute with each other and have 
eigenvalues ±1, a^^^ = 1, = I a^. Supercharges Qi of extended supersymmetric 
quantum mechanics form the algebra: 

{Q^,Qk} = 26ikH, {Qi,Qk} = {Q^,Qk} =0, i,k = 1,2, 

Qi = a^i\p + iVi+i{x)), Q, = a'i\p-iV^+i{x)), (2.1) 

where Vs{x) = Vi{x), = a± 1, = 1 ® a±. Hamiltonian and supercharges act on 
four-dimensional internal space and Hamiltonian is diagonal on vectors 'ip'^'^{x, E), where ai, 02 

are eigenvalues of (t\^\ cr^^^ ■ Supercharges Qi {Qi) act as lowering (raising) operators for in- 
dices o"!, (72. It is convenient to represent the Hamiltonian structure and the connection between 
wave functions in diagram form (see below). 

Obviously, due to commutativity of operators Qi and Qi with Hamiltonian, all %l)'^'^{x,E) are 
eigenfunctions of Hamiltonian with the same eigenvalue E. The only exception is the case, when 
the wave functions turn to under the action of generators of supersymmetry. 

Construction of isospectral Hamiltonians within A^ = 4 SUSY QM is based on the fact 
that four Hamiltonians are combined into the supermultiplet -ff^i- Nevertheless, it has to be 
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■tl:+(x,E) 



Q2 



Q2 



tp+{x,E) 



Qi 



i>+{x,E) 



noted that due to the symmetry under cxi -H- a2, i.e. H~_^^ = H"y-^_^ only three of them are 
nontrivial. Let's note, that such a relation in the case of higher-derivative SUSY [26] establishes 
the correspondence between a quasi-Hamiltonian and operators of the Schrodinger type, and is 
identical for any superpotentials. The procedure of construction of isospectral Hamiltonians, 
when ground state is removed from the initial Hamiltonian, is considered in |22j in detail. 
We will consider the construction of isospectral Hamiltonians by adding the states above the 
ground state of initial Hamiltonian. Similar procedure was already performed (e.g. in [23]), but 
the distinctive feature of the present research is to obtain of general results without specification 
of concrete form of the initial Hamiltonian. Let's consider the auxiliary equation: 



Hip(x) = eip{x). 
Let's take one of the Hamiltonians as the initial one 



H. 



0'2 



\{p - iaiV2{x)){p + i(JiV2{x)) + e = \{p - ia2Vi{x)){p + ia2Vi{x)) + e, 



(2.2) 



(2.3) 



where e is the so-called factorization energy. Hereafter the energy is measured from e. Strictly 
speaking, the supersymmetry relations with supercharges (|2.1|) expressed through superpoten- 
tial W{x) are satisfied for shifted by e Hamiltonian H — e. However, due to commutativity 
of supercharges with constant e, the relations between wave functions of H and H — e remain 
the same. When fixing operator H^^, the form of W{x) depends on the choice of factorization 
energy e, hence the Hamiltonians -ff^^ , H~^'^ , HZ"^ also have nontrivial dependence on e. 



When e < Eq (where Eq is the ground state energy of initial Hamiltonian), the auxiliary 
equation (12. 2p has two linear independent solutions (pi{x,e), i = 1,2, which are nonnegative 
and have the following asymptotics [15]: (pi{x) — t- +00 {(p2{x) — t- 0) under x — t- —00, and 
^Pi{x) — {(p2{x) — +00) under x — +00, i.e. the general solution has the form ip{x,£,c) = 
N{ipi(x,e) + Cip2{x,e)) with appropriately chosen constants (N is the normalization constant) 
and has no zeros. Thus, the function 'f{x, e, c) = ^-^ is finite and can be normalized at every 
concrete choice of e and c. Let's note that, for some values of e and c, ip{x, e, c) can have local 
extrema. In this case the natural choice of the initial Hamiltonian is or (which are 
identical due to the symmetry of H^^ under ai -H- (12). Then the superpotential has the form: 



W{x,e,X) 



-lln 

2 



1 + A 



dtif'^{t, e, c) 



Xi 



(2.4) 



with two new arbitrary parameters A, Xj, but one of them is inessential, because it gives an ad- 
ditional contribution to W{x). All the Hamiltonians forming the supermultiplet have nontrivial 
dependence on these parameters. 

To consider the connection between Hamiltonians from the supermultiplet, let's take as 
the initial one. Denoting the solution to 



m^p-{x,E) = Ei;-{x,E) 
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as ipj^{x^E) and using the first representation of the Hamiltonian H^'^ (the l.h.s. of (|2.3p ). we 
obtain the following relation between HZ, ipZ{x,E) and the initial expressions: 

_ ^ cf 
H_ = + —linp{x,e,c), 

1 W {i^^{x,Ei),ip{x,e,c)} 
w {x,Ei) = — z 7 

V'Z(x,£; = 0) = — - = (^(x,e,c). (2.5) 

if{x,e,c) 

The new state with E = Q and by definition normalized wave function (remind that energies 
are measured from e) appears in the Hamiltonian. For the discrete spectrum the normalization 
of the excited states wave functions is conserved. Using the second representation (the r.h.s. 
of (I23D) 

1 

2* 



= -{p-ia2V„A^)){p + ia2V„,{x))+e 



and identity Hj^ = the relation between H^, ip^{x,E) and initial Hamiltonian and wave 
functions can be obtained: 




<fix,e,c) 



^plix,E = 0) 



jV^V(x,£,c) 

{1 + Xf^dt^^t,e,c)y 



i;tix,Ei) = ^ ^ [ -^ + -^ln- ^]i^+ix,E,). (2.6) 

''+^ ' " y/2{E,-e) \dx dx {l + Xj^^dt^^{t,E,c))j''+^ ' ^ ' 

It is worth mentioning that the normalization of the wave function il^\{x, E = 0), as in the case 
of one-well potentials, can always be performed at any (^(x, e, c) and A by use of the following 
expression in the normalization condition: 

N-^^^{x,e,c) _ N-^ d 1 



(1 + AiV-2 J-^ rfi^2(i^ c))2 dx (1 + AiV-2 J-^ ^i^2(i^ £^ • 

From this relation it is easy to derive the relation between the normalization constants: N^'^ = 
(1 + A)A^~^. The normalization of ip^{x, Ei) is the same as for ip^{x, Ei) for any oi(p{x,e,c). The 
usage of superpotential (|2.4p with (p{x, e, c) corresponds to exact super symmetry, which leads 
to existence of zero- modes in HZ and H^. Existence of two zero- modes in super Hamiltonian 
of = 4 SUSY QM is caused by the fact that the Witten index theorem has to be modified 
when intertwining conditions are nonlinear, as discussed in detail in [25j . 

Let's consider the case when expression (|2.4p contains one of the particular solutions, e.g. 
ipi{x,e), instead of ip{x,£,c). If one of the particular solutions of second order differential 
equation is known, the second solution can be obtained from the relation 

P 1 

ip2{x,e) = ipi{x,e) / dt- 



>2(t,e)- 

Thus, the superpotential (12. 4p becomes 



W{x,e,\) = -\ln(l + X [ dt^^) 
2 V J-oo nit,i^)J 
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2 \(pi{x,e) + X(f2{x,£) J 2 \(p{x,e,X) 

It is easy to show, that the state with the energy = in the spectrum of HZ is absent (i.e., 
the wave function ipZixjE = 0) is nonnormalizable) hence the spontaneously broken N = 2 
supersymmetry exists. On the one hand, the zero energy state has the wave function 

V'|(x,£; = o) ^ 



ip{x,e,X) ' 

normahzable for certain values of A > 0, appears in the spectrum of H^. Here the exact N = 2 
supersymmetry takes place. Furthermore, it is known [29\ I30j. that the partial supersymmetry 
breaking is impossible in = 4 SUSY QM without central charges. This contradiction resolves 
with taking into consideration that the employment of a factorization energy e in the construction 
of isospectral Hamiltonians is the simplest way of incorporation of central charges in = 4 
SUSY QM. The similar situation occurs in consideration of form-invariant potentials [31]. More 
complete and consistent consideration of = 4 SUSY QM with central charges is given in [32], 
but this point is over the scope of the current paper. 



3 General properties of multi-well potentials 

We get started the consideration of properties of isospectral Hamiltanians, derived in previous 
section without appealing to the form of the initial Hamiltonian, from the following usefull 
expressions. Let yi and y2 be two linear independent solutions of a homogeneous second order 
differential equation. Then the following expression holds 

W{y„y2} 1 



/ A2yi{x)-Aiy2{x) \ _ / A2yi{xi)- Aiy2{xi) 
V Aiyi {x)+A2y2 (x) J \Aiyi (x, ) +^2?/2 (xi) 



Here W{yi,y2} = yiy'2~y'iy2 is the Wronskian, which for the second order differential equation, 
reduced to canonical form, as the Schrodinger equation, is independent on x and thus could be 
taken out of the integral. This expression is very useful for calculations of integrals in expressions 
First of all, it is natural to set Xj = — oo in (|2.4|) . because 



ip[x,e,c) 



(pi{x,e) + cip2ix,e) 

tends to in the limit and the function with asymptotic — ^ — ^ ~^ under x — t- — oo always 
exists, when we use a particular solution {pi{x,e). Therefore 

^~ / 9 = —r. ^KT777 7 [Mx, e, c - A -oo, e, c , 

= -(iTcWW^} " ' 

Aix,e,c) = '^^^^4^^f^. (3.1) 

<fi{x,e) +c(p2[x,e) 

Let's fix c = 1 in (j2.6p . This choice allows one to simplify the consideration, but nevertheless, 
reveals the fundamental features of Hamiltonians and wave functions of = 4 SUSY QM, 
using only general properties of the solutions to the auxiliary equation under e < Eq. Using 
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relations (|3.ip . it is easy to obtain the following modulo constant term expression for the super- 
potential: 



W{x,e,\) = --\n(l + \N-^^ T 

2 V J^oo{'fi{x,e) + Lp2{x,e)) J 



^ '^^^^'^^ + A(e, A)(^2(a;,e) 



2 V (pi{x,e) + ip2{x,£) 
. . ,^ ^ A(oo,£,l)-A-(A + l)A(-oo,£,l) 

^'^ A(oo,e,l) + A-(A + l)A(-oo,e,l)' ^ ^ 

Let's note that A(iboo,e,A), entering A(e, A), are determined by the asymptotics of solutions 
to the auxiliary equation. Due to this fact, and since for HZ the potential is determined by 
the symmetric combination (pi[x,£) + ip2{x,£), while in the case of - by the asymmetric 
combination ipi{x^e) + h.{e,X)ip2{x,£), we get: 

_ _ (fi 
H_ = H^- -^ln((^i(j;,e) + 932(2;, e)), 

Ht = H^ - ^ln{ipi{x,e) + A{£,X)^2{x,£)) . (3.3) 

In some sense, the potentials UZ{x,e) and U^{x,e, X) are form-invariant [H], i.e. potentials 
and wave functions transform to each other by changing of parameters and their spectra are 
identical and this holds independently on the choice of the initial Hamiltonian. As it will be 
shown in the next section, if the potential in HZ is a multi-well symmetrical potential, then in H^ 
it should be asymmetrical. Moreover, varying e, as well as A, the form of U^{x, e, A) can change. 
Relations ()3.ip and ()3.2p are also useful to derive the exact form of the wave functions ipZix, E) 
and ip^{x,E). Thus, the expression for the wave functions 'il:^{x,E) comes from the similar 
expression for ^Z(x, i?) by the substitution 931(2;, e) + 932(2;, e) — )■ {(pi{x,e) + A{e, X)(p2{x,e)). In 
particular, the normalization constant for ip^{x, E = 0) can be obtained from the corresponding 
expression for ipZ {x,E = 0). 

To demonstrate the possibility of partial supersymmetry breaking in = 4 SUSY QM we 
have to show that the wave function il)^{x,E = 0) ~ ^^J^ is normalizable. The value of the 
normalization constant can be derived from (j3.ip under c = A > 0. We will return in what 
follows to the calculation of this constant for the concrete form of the initial Hamiltonian. In 
its turn, the potential in H^ has the same form as the corresponding potentials in the case of 
exact supersymmetry under replacement A(e, A) — t- A. 

As it follows from (j2.6p and (j3.3p . the parametric dependence on A is only present in H^ 
and ip^{x,E). It should be noted that the value of A is not to be fixed by the normalization 
condition for the wave function. Existence of a parametric freedom is common in building the 
isospectral Hamiltonians, which is based on different versions of the inverse scattering prob- 
lem [331 ESI [36] • It is caused by the ambiguity in reconstructing quantum- mechanical potentials 
from the spectral data. We will show latter that similar situation exists also in stochastic mo- 
dels. This is unexpectedly enough, since the potentials of stochastic models have direct physical 
meaning in contrast to quantum mechanical potentials. 

4 Probability densities functions 

The FP equation is equivalent to the Langevin equation, however its application in physics 
is more wide, since it is formulated in more appropriate language of the probability densities 
P^{x,t; xo,to). According to [9l|T2] the FP equation has the form: 

d D d 

—P^{x,t;xo,to) = —-^P^ix,t;xo,to) T -^V{x)P^{x,t;xo,to), 
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P^ix,t;xQ,to) = {6ix-Xo)) , U±{x) = ± dzV{z), 

Jo 

where U±{x) is the potential entering the Langevin equation. Let's set to = 0. The Fokker- 
Planck equation describes the stochastic dynamics of particles in potentials U+{x) and C/_(x) = 
—U+{x). Substituting 

P±(x,t;xo,0) =exp|--i[C/±(x) -i7±(xo)]|K±(x,t) 
the FP equation transforms into the imaginary time Schrodinger equation: 



-D^K^{x,t) = {-^^ + \ [V\x) ± D V'{x)] } K^{x,t) 
K±{x,t) = 



exp 



D J 



xo), Hi = -^^ + ^-[V\x)±DV'{x)], (4.1) 



in which the diffusion constant D can be treated as the "Planck constant", while H± has the 
form of = 2 SUSY QM Hamiltonian. Equations (|4.ip are basic to apply of the eigenfunctions 
expansion method for construct exactly solvable stochastic models by use of the results of cor- 
responding quantum mechanical problems. Details of this method in the framework of = 2 
SUSY QM can be found in PITS]. Utihzation of the extended N = 4 SUSY QM formalism gives 
additional possibilities to construct new exactly-solvable models of stochastic processes, since 
the Hamiltonian H^^ of N = 4 SUSY QM includes four isospectral Hamiltonians. It allows to 
obtain new K^^{x,t) and U"'^{x), hence P"^{x,t;xQ). 

To obtain the expressions for probability densities P^^{x,t;xQ) we should to know not only 
the wave functions and the spectrum of H^^, but those of the corresponding U^^{x). It is easy 
to see that 

C/+(x,c = 1) = -Dln(v9i(x,e) + (^2(x,e)), UZ{x,c = I) = -U:^{x,c = I). (4.2) 

Further consideration is based on account of the symmetry of = 4 SUSY QM H^^ under 
o"! 0(72, that leads to the relation: 

H^(x,p) = ^QS^'Q^^ ^ ^Q^^^O^^^ = H+{x,p). (4.3) 

First equality indicates the existence of the expression of H^'^^ in terms of Qi (Qi), while 
the second equality implies that H^~^ is expressed in terms of Q2 {Q2)- At the same time, the 
supercharges entering i^^^-* and are substantially different: 

H+{x,p) = ^Q(-^)q(+) = \[p' + {vl'-\x)f-DVI'-^'{x)], 



Vl^\x) = D^ln 
ax 



<fix,e,c) 



1 + A f!^, dx' [(^(x',e,c) 



-D^\u\^i{x,e)+ K{e,\)^2{x,e)\ . (4.4) 



According to (14. 3p . quantum H_ and have the same spectrum and the wave functions, 
thus the corresponding K{x, t) are identical. But the corresponding stochastic models are de- 
scribed by essentially different potentials, therefore their P(x,t;xo) are different. Moreover, the 
potential U'^{x,£,A) = —U^{x^e,K) = —Dhi\Lpi{x,£) + A(e,A) ip2{'^,£)\ has non-trivial para- 
metric dependence on A, which implies the existence of the family of stochastic models with the 
same time dependence of the probability density, but substantially different coordinates depen- 
dence. Surprisingly, the parametric freedom allows for a modification of the shape of potential. 
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Recall that physical quantities, e.g. time of passing the potential maximum, considerably depend 
on the local modification of the potential form [371 ESI 132] • 

It is easy to obtain the corresponding P^^{x,t;xo) using the expressions for wave functions 
ip^l{x,En) and potentials U^^{x) from ([23]), (pIB]) and 1^^ . (g^]). It is important to note, 
that transitional probability densities {x,t; xq) and P^{x,t;xo) have also the parametric 
dependence on A. 

Denoting ip{x,e,c = 1) = (pi{x,e) + (p2{x,e), the expression for calculation of distribution 
function P~(x,t;xo) takes the form: 



PI{x,t-xo) = ^ , + (^^^ ) 5^e-i^(^-^)^I(x,i?„) ^PZ{xo,E,,). (4.5) 

{^{x,e,l)y \ip{x,s,l) J 

For simplicity let's assume that the spectrum of the initial Hamiltonian has only the discrete 
states. When the continuous spectrum exists the summation should be replaced with integra- 
tion over the corresponding density of states. Substituting the expressions for wave functions 
tpZ{x,En) into (jl3D we obtain: 



TV 



1 



P_ {x, t; xo) = — —2 + - {(p{xo,e, l)y 



(En-e) dxdxo yip{x,e,l) J y(p{xo,e,l) J' 

The expression for P^{x, t; xq) can be obtained from equation (|4.5p by substituting ip{x, e, 1) — 
ip{x,e,A{X,e)), N'"^ NJ[^. Both of Pl{x,t;xo) and P+( x,t;xo) have the equilibrium distri- 
butions at t — >• oo, which coincide with square of the zero modes wave functions in Hamilto- 
nians HZ and H^. Moreover, it is easy to see that terms from the excited states of HZ and H^ 
do not contribute to the normalization condition of these density functions, i.e. the normaliza- 
tion of these functions is guaranteed by their equilibrium values. Dependence on parameter A 
in P^ {x,t; Xq) is not eliminated by the normalization condition. 

The distribution P^ (x,t; xq) {P^{x,t]Xo)) describes the stochastic dynamics of a particle in 
the metastable state of the potential, inverted to U^{x,e) (— C/^(x, e, A(e, A)). Thus, it has not 
the equilibrium limit and is determined by the expression: 

(x 1) ^ E - t 

P+{x,t;xo) = '^^ ^\ y^e-^-^^i;^{x,En)i^+{xo,En). (4.7) 



99(0:0, e, 1) 



n=0 



Similarly to the previous case Xq) can be obtained from P^{x,t;xo) by substitution 

ip{x,e,l) — )• ip{x,£,A{\,e)). The parametric dependence on A can be eliminated by the nor- 
malization condition for P^ {x,t] xq), but as it will be demonstrated below, choosing the initial 
Hamiltonian of harmonic oscillator, the parametric arbitrariness in P^ {x,t; xq) remains. 

To summarize, we should note that in the framework of = 4 SUSY QM four models of 
stochastic dynamics emerge. Two of them correspond to the motion of particle in double-well 
potentials (symmetric and asymmetric) under the action of the Gaussian noise, the other ones 
describe the dynamics of the particle in the metastable state, when iV = 4 supersymmetry 
is exact. For the case of the partial supersymmetry breaking PZ {x,t; xq), P^ {x,t; xq) and 
P^{x, t; Xq) describe stochastic processes in the metastable state, while P^{x, t; xq) - in bistable. 
Some of them possess the parametric freedom, that allows to change the parameters of external 
field. 
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5 Isospectral Hamiltonians with almost equidistant spectrum 

To construct the explicit expressions for potentials and wave functions we choose the initial 
Hamiltonian with the harmonic oscillator (HO) potential. Let's consider the solution to the 
auxiliary equation (or e < Eq = ^, {h = m = 1): 

£ . 



Introducing dimensionless variables ^ = we obtain the equation for where e 

Now we consider two different cases. 

a) Exact AT = 4 SUSY. This equation has two linear independent solutions: parabolic 
cylinder functions Dy{\f25^)^ Dy{—^/2^). According to terminology in the above, we denote 
ipi{i,e) = Dy{V2C), ip2{i,e) = D^^-V^i) and the Wronskian becomes W{ipi,'f2] = 00], 
with gamma-function r(— i^). Following the procedure from the previous section, the general 
solution to the auxiliary equation is chosen to be: 

ip{i,e,l) = D,{V2ei+D,{-V2ei. (5.1) 

As one can see from (jS.ip . (p{x,e, 1) is an even function of ^. To obtain the exact form of the 
superpotential, the integral, entering the definition of the W{x, e, A) as well as the normalization 
constant N^'^ have to be calculated from p.ip . (j3.2p . Due to the symmetry of (/?(^, e), expression 
of the integral simplifies and takes the form: 



1 + XN^^ 



dt ^ XN- 



,{<fl+(p2f 2W{ipi,ip2} 

Hence, the superpotential is: 



991 (g,£) - (y?2(^,e) 



+ A(+oo,e, 1) 



\ 'Pi + ^2 J 

^-2 _ W{(pi,ip2} _ 2/? 



A(+oo,e,l) T{-iy) 

Using the asymptotic value of the parabolic cylinder function, we have A(+oo,e, 1) = — 1 and 
the expression for the superpotential simplifies: 

From (15. 3p it follows, that the value of A is restricted to A > —1. Using (12. 5p for HZ and 
ipZ{x,E), we can derive the form of the Hamiltonian and the corresponding wave functions 
(here 'i/'+(^> -E'j) are the wave functions of the HO): 

HZ=H~-^\n (l).(V20 + D,{-V2i) 

"^-^^^^'^ = vmrf) — ^i^) — ' 
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Figure 1. Shape of the potential U_ e). Figure 2. Potential ^/^(C, £ = 0.47, A). 



Analysis of (j5.4p reveals that there exist several local minima only for < e < ^. With values 
of e to be close to the right boundary, the third local minimum appears (Fig. [T]) and depth of 
the outside minima increases. It should be noted that in terms of the dimensionless variable ^ 
the only way to vary the form of the potential is by varying e and A. In the case of natural 
units, additionally, the form of the potential (in particular, positions of the local minima) can 
be changed by variation of u. 

Connection between H^, ip^{^,En) and H^, i/^^(^,£'„) can be viewed in the same manner. 
Expressions for and 'iJj^{^,En) can be obtained from (15. 4p by substituting (^(^,e, 1) — )• 



HX = H^-^ In {D,{V2i) + (1 + \)D,{-V20) , 

1 t^{^;(^,i?0,<^(e,e-,A + i)} 



yj2{E, - e) (^(^,e-,A + l) 



{D,{V20 + (A + l)D,i-V20) V>{x, s, A + 1) ' 
2 _ 2(A + 1)V¥ 

" r(-z.) ■ 

From these relations it follows the restriction — 1 < A. As it can be seen from Fig. [2l the 
potential (Ci A) possesses the well indicated asymmetry, which increases under A — ?• — 1. The 
value e = 0.47 corresponds to the region, where UXiC^^) A) has three local minima. It should be 
noted, that the depth of the central minima also increases with A — )• — 1. The obtained potentials 
and the corresponding wave functions may be applied for the resonant tunneling phenomenon [3] 
studies. 

b) Partial N = 4 SUSY breaking. The utilization of the particular solution of auxiliary 
equation for derivation of superpotential (j5.2p realizes the situation of partial supersymmetry 
breaking. As it was mentioned above, the spectrum of HZ does not contain states with E = 0, 
because their wave function is nonnormalizable. The shape of UZ is presented in Fig. [3l 
At the same time, zero state appears in the spectrum of with the wave function: 

+ . X NZ^ 



(^(e,e,A) (^i(e,e) + A(/P2(^,e)) 
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Figure 3. Potential UZ{£„e) (0 < £ < 0.5). 



The normalization constant ^ is calculated by use of (j5.4p and has the form ^ = j^—^- 
Thus, it is restricted to A > 0. The Hamiltonian becomes = — ^ In ((^(^, e, A). The 
potential shape {e = 0.47) is presented in Fig. [2] with A > 0. 

When e = -i (D_i(V20 = e" (1 - ^>(e)), with $(0 to be the error function) , potentials 



U_ (x, e) correspond to that of [H] (the only difference is in additional constant term) 



l\ 4 

+ 



^(1-^(0) 



^(1-^(0) 



(5.5) 



This potential is single- well, as well as corresponding ?7^(^, e = — ^,A), because existence of 
several local minima is possible only for < e < ^. Using the form-invariance property of U^: 



-^ln((l + A)-(l-A)$(0). 



(5.6) 



The spectrum of the Hamiltonian with ()5.6p contains, in contrast to (|5.5p . the state with E = 
and the wave function of the form: 



((1 + A)-(1-A)cl>(0)- 



6 Generalization of the Ornstein— Uhlenbeck process 

Let us give an example of the construction of new stochastic models. We choose the harmonic 
oscillator (HO) Hamiltonian as the initial one. Then, we consider the auxiliary equation with 
e < £"0 = (recall, that D plays the role of the "Planck constant" in the FP equation and 
e = T^j). Here we consider two cases. 

a) Exact N = 4 SUSY. For many problems of stochastic dynamics it is reasonable to 
use symmetric double- well potentials, so we consider first the case c = 1. As it can be seen 
from (j5.ip . ip{^,e,l) is even function, and the normalization constant N~'^ can be calculated 
from (j5.2p . One of the main characteristics of the stochastic process is the potential entering the 
Langevin equation. As it follows from ()4.2p . UZ{x) = —Ul^{x) = Dlinp{x,£,l) is symmetric, 
while U^{x, A) = —U^{x, A) = D In ip{x, e, A + 1) is asymmetric, that depends on the value of A. 
The existence of the parametric freedom in U^{x,A) is quite surprising: contrary to quantum 
mechanics, where potentials reconstructed from the spectral data indeed have a parametric 
freedom, but give the same observables, in stochastic mechanics local properties of the potential 
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U 





Figure 4. Potentials U^l- a) exact TV = 4 SUSY: UX{£„\) {oj = D = 
N = A SUSY: U^{^) - lilac, UZ{£,) - blue, U+{t A) - dashed, C/+(e, A) 



1, £ = 0.47); b) Partial broken 
- red (e = 0.47, A = 1). 



substantially influence on many characteristics, for example, on the rates of the potential barrier 
crossing [371 [381 [39] . 

U'^{x,K) = —U^{x,K) = — Z) In e) + (A + 1)992(3^, obtained within the framework 

of proposed model, is presented in Fig. dj Changes of the U^{£^,\) shape due to varying A 
become more and more sharp under A — t- —1. The analysis of the expressions of U^{^, A) shows, 
that several local minima are possible only when < e < ^. Variation of e inside this range 
changes the height of the barrier, which significantly increases with e — )• ^. Moreover, the shape 
of the potential can be changed (especially the minima locations) by varying w, when switching 
to natural variables. To summarize, the existence of (w, e. A) significantly change the shape of 
the C/+(x,A). 

Substituting the HO wave functions to ()4.7p and using the Mehler formula |40j . we get 



(6.1) 



where z = e As it was noted in the above {i^z;^^) does not have the equilibrium value 
and could be normalized (the proof of this claim is given in Appendix), i.e. 



/+00 
diP^{i,z-i^) = l. 
'OO 



The probability density P_ {£,,z;^q) is characterized by the existence of the equilibrium value, 
which is determined by the zero energy wave function of HZ, i.e. 

PI {x,t ^ 00; xq) ^ _ iV"^ - 



ipHx,e,l)' r(-i.)- 
After some calculations (see Appendix) we obtain the expression for PZ {S,, z; S,o)- 

m Co, r) = (-^ [F{^, ^0, r)] ^{^0, e, 1) - F{C, Co,r)-^ [ifi^o, e, 1) 



F{i,io,T) = exp 



2 1 - 
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The probability density Pl^(^,z;^o) can be obtained from P^{^, z;^o) ()6.ip by substituting 
(p{x,e,c = 1) —7- (p{x,e, X + 1). Note that the normahzation condition does not ehminate 
the A-freedom (see Appendix). As it has been noted in previous section, the expression for 
P^{^, z;£,o) is obtained from PZ{C,z;^o) with replacing ip{x,e,c = 1) — > ip{x,£,X + 1), A*""^ — )• 

- r(-.) • 

The dependencies of P^{^,z]^o) on spatial and time variables, presented in Fig. [5l demon- 
strate sharp modifications (with fixed e, uj, D) at different values of A, especially with A — t- — 1, 
when the considerable assymmetry of P^{^,z;(^q) as well as decreasing the time of passing to 
bimodality are observed. Moreover, modifications of the shape of potential U^{^, A), as well as 
-P+ (^, z;^o), can take place with variation of oo, when we pass from dimensionless variables ^ to 
physical x. This leads to the shift of the local minima and change the height of the barriers. 
When e — )• |, the barrier between local minima significantly increases, that allows to study the 
transition time from one local minimum to another as a function of 

b) Partial N = 4 SUSY breaking. Let's briefly discuss partial A^ = 4 SUSY breaking 
and its consequences to stochastic dynamics. In this case three models for the description of the 
stochastic dynamics of particles in metastable states {U^{^, e), UZ{^, e), U^iC^ A)) and one in 
the bistable state (?7^(^,e, A) ) arise in the considered approach. The corresponding potentials 
have the form: 

UZ{i,e) = -U^{i,e) = D\n\D,{V2i% 

UX{i,e,\) = -U^{i,e,X) = D\n\D^{V20 + \D^{-V2i)\. 

The range of A variation in J7^(^,e, A) and C/j^(,^,e, A) is restricted to be A > 0. Fig. 
presents the dependencies of the potentials on spatial variables at fixed e and A. Stochastic 
dynamics in mentioned metastable states satisfies to different boundary conditions: reflective 
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boundary at left (right) and absorbing at right (left) for UZiC^S) {U^{(,,e)) and absorbing 
boundaries at both sides for U^{^,£,X). The corresponding probability densities can be calcu- 
lated with taking into account previously obtained expressions with appropriate modifications. 

7 Conclusions 

In the paper we consider the construction and the study, within the framework of = 4 
SUSY QM, of the general properties of Hamiltonians with multi-well potentials. These results 
have been used to obtain exactly solvable models of stochastic dynamics. The special attention 
was made to studies of features of these Hamiltonians without concretization of the form of 
initial Hamiltonian. Relations for certain type of integrals, containing the fundamental solutions 
to the Schrodinger type equations, allows one to show the form-invariance of the isospectral 
Hamiltonians with multi-well potentials, obtained within of = 4 SUSY QM. In other words, 
having the identical spectra of HZ and H^, the corresponding potentials can be obtained from 
each other through replacing parameters of the problem. Moreover, this relation allows to 
analytically calculate the normalization constants of zero modes, using only asymptotic values 
of fundamental solutions. The construction of isospectral Hamiltonians within = 4 SUSY QM 
also implies the partial supersymmetry breaking. Taking the model of harmonic oscillator as an 
example we derive the exact form of the isospectral Hamiltonians with multi-well potentials and 
the corresponding wave functions. The existence of a parameteric freedom gives a possibility 
to vary the shape of the potential in the wide range. This becomes important in the research 
of different phenomena, such as tunneling processes, which are sensitive to the structure of the 
multi-well potentials. 

The obtained results have been used for the construction of exactly-solvable stochastic models. 
Obtained potentials, entering the Langevin equations, and probability functions are characteri- 
zed by the parametric dependence, which allows sufficiently modify their shape. A parametric 
freedom is typical for quantum mechanical isospectral Hamiltonians, where different versions of 
the inverse scattering problem are used [341 [35| l36] . Such a freedom arise in view of an ambiguity 
reconstructing the potential from the spectral data. On the other hand, the parametric free- 
dom in stochastic models is quite surprising, since their potentials, in contrast to the quantum 
mechanics, have direct physical meaning. Many characteristics of stochastic dynamics, such as 
potential peak passage times and the metastable state lifetime [37 1 138 ^ 139]. substantially depend 
on the shape of the potential. Therefore, it is interesting to investigate the characteristics of 
stochastic processes in multi-well potentials, such as the Kramers problem, the stochastic reso- 
nance etc., and their dependence on a modification of the shape of the potential. We note that 
our results admit a generalization to non-Markovian processes, by use of the approach of |42j . 
Furthermore, the obtained results can be generalized to the stochastic models of polymer dy- 
namics [l3]. We have also considered non-trivial consequences of the partial supersymmetry 
breaking in = 4 SUSY QM to the description of stochastic dynamics. In this case three 
probability density functions with the same time dependence (i.e. with the identical spectrum 
of the Fokker-Planck operator) describe stochastic dynamics of particles in metastable states, 
and one function corresponds to the dynamics in the bistable state. 

A Appendix 

Let's consider details of obtaining the expression of Pi z; ^q) and z; (,o). According 

to (14. 6p the distribution function PZ{(^, z;(,o) has the form: 



P_ {x,t;xo) 



2 + 2 (v'(^o,e, 1))^ 



1)) 
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For the generalized Ornstein-Uhlenbeck process, substituting (jS.ip and the wave function of 
harmonic oscillator il)^{^,En) = (;^)^^^ (jA.ip we get the expression: 

1 / g; \i/2 d 1 ^ 

X Z)„(V2e) (^^Z)„(V2^o)^(eo,e, 1) - A.(\^eo)^9'(eo,e, 1)) , 
where ^ = y^x. Using the Mehler formula |40j : 

2^ -,Dn{x)D,,{y)z = -== exp ^— ^JZ^ = ^(^' ^) 

n=0 ■ V i ^ 1. J 

and the relation 



V \ -Z5„(x)Z?„(y) = / drT"-(5+-")F(x,y,r) 

^ n! (n + i - e) Jo 



we obtain the expression for PZ {C, z; ^o) (!6.2p . It should be noted that presence of spatial deriva- 
tives eliminate the contribution from excited states of the Hamiltonian HZ to the normalization 
condition for PZ {S,, z; £^q). 

Let's consider the question of normalizability of P^ {x,t; xq) and P^ {x,t; xq) as well as the 
possibility of elimination of the A-freedom in PI^{x,t;xQ) by the normalization condition. As 
it was noted in the above, P^ {x,t] xq) is obtained from P^(x,t;xo) by replacing ip[x,e,l) — >• 
if{x, e, A + 1) and has a form 



^- ^' Wd) e-^VMio.e, A + 1) ^/^^ 
The normalization condition for the distribution function: 

dxP+{x,t;xo) = J- / dCP+{^,z;^o) = l- 

-oo V W 

Substituting the integral representation of the parabolic cylinder function 



z2 

Dp{z) = y J e-'^^'-x-P-^dx, Rep < 0, 



results in 

^2 



dC e lT^e-« /2 L'^(±^/2e) = z\\ - 2^)1/2 e-«o/2^^(±^^p), (A.2) 

-oo 

Substituting (|A.2p to the normalization condition for Pj~(x,t;a;o) and P_p(x,t;xo) it is easy 
to verify, that they are equal to one, thus the probability densities, which correspond to the 
stochastic dynamics in a metastable state, are normalizable and the normalization condition 
does not eliminate the A-freedom in t; xq). 
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